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ABSTRACT 

The dust shell around the evolved star HD 179821 has been detected in scattered 
light in near-IR imaging polarimetry observations. Here, we subtract the contribution 
of the unpolarized stellar light to obtain an intrinsic linear polarization of between 30 
and 40 per cent in the shell which seems to increase with radial offset from the star. 
The J and K band data are modelled using a scattering code to determine the shell 
parameters and dust properties. We find that the observations are well described by 
a spherically symmetric distribution of dust with an density law, indicating that 
when mass-loss was occurring, the mass-loss rate was constant. The models predict 
that the detached nature of a spherically symmetric, optically thin dust shell, with a 
distinct inner boundary, will only be apparent in polarized flux. This is in accordance 
with the observations of this and other optically thin circumstellar shells, such as 
IRAS 17436-1-5003. By fitting the shell brightness we derive an optical depth to the 
star which is consistent with V band observations and which, assuming a distance of 
6 kpc, gives an inner shell radius of r-jn — 1.44 x 10^^ m, a dust number density of 



2.70 X 10 m at Tin and a dust mass of Md — 0.08 Mq. We have explored 



axisymmetric shell models but conclude that any deviations from spherical symmetry 
in the shell must be slight, with an equator-to-pole density contrast of less than 2:1. 
We have not been able to simultaneously fit the high linear polarizations and the small 
{E{J — K) = —0.3) colour excess of the shell and we attribute this to the unusual 
scattering properties of the dust. We suggest that the dust grains around HD 179821 
are either highly elongated or consist of aggregates of smaller particles. 
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1 INTRODUCTION 

The evolved star HD 179821 (IRAS 19114+0002) is sur- 
rounded by a dusty circumstellar envelope which is de- 
tectable in scattered light out to a radius of at least 9 arc- 
sec (Kastner & Weintraub 1995; Ueta et al. 2000). Mid-IR 
imaging (Hawkins et al. 1995) shows that the warm dust is 
distributed in a shell around the star with an inner bound- 
ary at a radius of 1.75 arcsec (Jura & Werner 1999). There 
appears to be little or no dust inside this inner boundary, 
suggesting that the mass-loss terminated some time ago. 
The detached nature of the shell was also indicated by ear- 
lier ground-based and IRAS photometry which showed the 
double-peaked profile characteristic of emission from a hot 
photosphere and warm dust shell (Hrivnak, Kwok & Volk 
1989; hereafter HKV89). 

It is still not clear whether HD 179821 is a low or inter- 
mediate mass post- AGB star or a more massive supergiant 
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(see Reddy & Hrivnak 1999 for a review of the arguments, 
also Thevenin et al. 2000). This fundamental uncertainty is 
reflected in the wide range of estimates of the distance to 
the star (between 1.5 and 6 kpc) and the quantities depen- 
dent upon it, such as the physical size of the shell and the 
mass of dust contained. If the star is in the post- AGB phase 
then the object is a proto-planetary nebula (PPN) and sim- 
ilar to other optically thin shell-like PPNe such as IRAS 
17436+5003 (Ueta et al. 2000; Gledhill et al. 2001) 

Near-IR imaging polarimetry of HD 179821 (Gledhill et 
al. 2001) allowed the overlying stellar point spread function 
(PSF) to be separated from the scattered light to provide an 
image of the dust shell in polarized flux in the inner 4 arcsec 
around the star. These images show a spatially resolved ring 
of scattered light encircling the star, with an inner radius of 
between 1.1 and 1.6 arcsec (the star is not at the centre 
of the shell), indicating that both the near-IR and mid-IR 
emission arise from the same body of dust illuminated and 
heated by the star. By subtracting the stellar PSF, an image 
of the shell surface brightness in the near-IR, and hence the 
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Figure 1. A J band image in polarized flux of HD 179821 (IRAS 
19114+0002). The data are from Gledhill et al. (2001). A circle 
of radius 1.9 arcsec, passing through the regions of peak flux, 
illustrates the symmetry of the shell. The centre of the circle is 
denoted by +, which is offset from the stellar centroid, denoted 
by *. 



intrinsic shell polarization, can be formed. In this paper we 
use these quantities, along with the distribution of polarized 
flux, to fit scattering models to the observations in order 
to deduce the nebula geometry, dust distribution and the 
properties of the dust grains. 
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2 OBSERVATIONAL CONSTRAINTS 

2.1 Imaging polarimetry observations 

A J band image of HD 179821 in polarized fiux, obtained at 
the 3.8-m UK Infrared Telescope (UKIRT), is shown in Fig. 
1. The polarimetric observations were made with the IR- 
CAM infrared detector and the IRPOL polarization module 
and are described by Gledhill et al. (2001). A 2x warm mag- 
nifier gave a pixel scale of 0.143 arcsec. Fig. 1 clearly shows 
the scattered Ught to lie in a shell or ring surrounding the 
star. A circle of radius 1.9 arcsec is drawn on the image to 
indicate the peak intensity (in polarized fiux) in the shell. 
Although there is evidence for structure in the scattered 
light, to a good approximation the shell is symmetric, with 
equal amounts of fiux lying inside and outside the circle. The 
centre of symmetry (+), is offset from the stellar centroid 
(*) which lies 0.3 arcsec to the N, indicating that the star 
does not lie at the centre of the shell. Jura & Werner (1999) 
noted this offset in their 11.7 /im images which show that 
the thermal emission is also distributed in a shell around the 



Figure 2. Azimuthally averaged radial profiles through the shell 
in the J (upper) and K (lower) wavebands. Polarized flux mea- 
surements are shown as crosses and are plotted on an arbitrary 
flux scale. Per cent polarization is plotted as triangles according 
to the scale on the right. 

star, as first observed by Hawkins et al. (1995). The bright 
shell seen in polarized fiux is surrounded by a fainter halo 
of scattered light, which can be seen in our images out to 
a radius of approximately 5 arcsec from the star. In addi- 
tion. Fig. 1 shows that polarized flux is recorded inside the 
bright shell. The brightness of HD 179821 in the near-IR (K 
— 4.59 ±0.05, Kastner & Weintraub 1995) means that align- 
ment and subtraction residuals dominate the central region 
of our images, so that we have not been able to measure 
polarization closer than 1 arcsec to the star. 

We take advantage of the radial symmetry to form az- 
imuthally averaged radial profiles of polarized fiux and per 
cent polarization, shown in Fig. 2. The profiles are centred 
on the + symbol in Fig. 1. Since the observations were made 
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through varying degrees of cloud cover, it has not been pos- 
sible to calibrate the flux scale. The peak in polarized flux 
occurs at a radius of 1.9 arcsec at J and 2.0 arcsec 
Given that the pixel scale is 0.14 arcsec and that the 
conditions at the time were variable, the difference 
significant. Per cent polarization peaks at 18.5 per > 
J and 18.0 per cent at K, both at a radius of 2.7 
In fact, the radial profiles in both polarized flux a 
cent polarization are the same, within errors, at J < 
This suggests that the scattering properties of the d 
remarkably wavelength independent. 

2.2 PSF subtraction 

The polarizations plotted in Fig. 2 were obtained 
viding the polarized flux (Fig. 1) by the total flux 
However, this includes a contribution from the unpc 
point spread function (PSF) of the star, whose wii 
tend to cover the whole of our image. These polar! 
are, therefore, lower limits on the intrinsic polarizatioi 
shell, Pint. To estimate the intrinsic polarization, w 
subtract the PSF of the star to obtain an image of tl 
in total flux so that 



'psf 



where, /, Q, U are the Stokes intensities in the unsub 
data and 7psf is the subtracted PSF flux. Note that i 
larization angles (and hence the centro-symmetric sea 
pattern) depend only on Q and U and are not affec 
the PSF subtraction. 

To perform the PSF subtraction we used scali 
aligned images of nearby bright stars: HD 178744 
V=6.3) for the J band and HD 179626 (F7V, V=c 
K band. Both stars were observed in the same poli 
ric mode as the target for a total of 60 sec each, re 
in a signal-to-rioise ratio (SNR) in the wings of the 
greater than 2 out to a radius of 4 arcsec. The scali 
based on the peait pixel values of the stars. The PSF si 
tion can only be approximate; seeing conditions at tl 
were variable and no better than ^ 1 arcsec FWHM. ] 
in all of our subtractions we found that the central pi 
surrounded by a negative 'hole' indicating over-subtrauuiuu. 
This suggests that either (i) the alignment between the tar- 
get and PSF star could not be achieved to sufflcient accuracy 
to correct the central bright regions or (ii) the shape of the 
PSF star was slightly different in the central peaJj (broader) 
compared with the target. Both could be due to slight vari- 
ations in seeing. A reduction in the subtracted flux of 10 per 
cent was sufficient to turn the negative hole into a positive 
pedestal, so that the central flux is then under-subtracted. 
Subtraction of the PSF approximately doubles the polariza- 
tion in the shell in both wavebands to between 30 and 40 
per cent, so an error in the PSF subtraction of ±10 per cent 
represents an error of ±5 per cent in the intrinsic polariza- 
tion. The fact that the polarization is observed to double 
indicates that although the SNR in the wings of the stan- 
dard star is low compared to the target, there is sufficient 
signal to correct for the stellar PSF. 

After subtraction of the PSF, radial profiles were formed 
of the total flux in the shell and of the intrinsic per cent po- 
larization and these are shown in Fig. 3. The polarization 
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Figure 3. Azimuthally averaged radial profiles through the shell 
in the J (upper) and K (lower) wavebands, after subtraetion of 
the stellar PSF. Total flux measurements are shown as eircles and 
are plotted on an arbitrary scale. Per eent polarization is plotted 
as triangles according to the scale on the right. At radii less than 
1 arcsec, the total flux becomes negative indicating that the PSF 
has been over-subtracted in the central regions and consequently 
the turn over in total flux at 1.8 arcsec radius may be an artifact. 
The dashed curve on the upper plot is V band HST data (see 
text). 



now rises with increasing offset from the star, suggesting 
that the depolarizing effect of the stellar PSF increases with 
radius. As with the polarized flux, the profiles of total flux 
and per cent polarization show little variation between J 
and K. The error bars shown in Fig. 3 are derived from 
the spread in pixel values during the azimutlial averaging 
process and do not include any systematic errors that may 
arise from the PSF subtraction (above). In Figure 3 we also 
plot the profile from the V band HST image of Ueta et al. 
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which we have obtained by smoothing their image 
to the 0.9 arcsec FWHM seeing of our UKIRT data, forming 
the azimuthaUy averaged profile and then scaling it to our 
intensity profile at a radial offset of 1.9 arcsec. Although the 
HST image has not been corrected for the WFPC2 PSF, 
its contribution to the shell brightness at several arcsec off- 
set from the star will be much less than in our ground- 
based data. The similarity of the HST and PSF-subtracted 
UKIRT profiles gives us confidence that our PSF subtrac- 
tion is fairly sound. Indeed, both profiles approximate an 
surface brightness dependence, as expected for optically thin 
scattering in a shell with an density fall-off. At a radial 
offset of less that 1.5 arcsec, the PSF-corrected fiux shown 
in Fig. 3 starts to decrease, and at radii of less than 1.0 arc- 
sec, it becomes negative. This indicates an over-subtraction 
in the central regions. There is no indication from the HST 
data of a decrease in shell brightness in the central regions, 
and so this may be an artifact of the PSF subtraction. It 
is also possible that too much flux has been subtracted at 
radii > 1.5 arcsec and that the intrinsic polarization is lower 
than we have estimated. To try to further corroborate the 
PSF subtraction, we also used a separate J band observa- 
tion of HD 179821 and PSF standard, taken in June 1999 
(a year after the polarimetry observations) to form a PSF- 
subtracted image of the shell. We then divided this image 
into our (May 1998) polarized flux image and find that the 
shell polarization is the same, within errors, to that shown 
in Fig. 3. 



2.3 Shell properties 

2.3.1 Relative brightness and colour 

A measure of the scattering cross-section, Csca, of the dust in 
the shell is provided by the ratio of the shell surface bright- 
ness to the brightness of the star. This ratio, which we call 5", 
is intrinsic to the shell and independent of interstellar extinc- 
tion. The variation of 5* with wavelength provides a strong 
diagnostic of the dust model, in particular the size distribu- 
tion. To quantify the surface brightness of the shell, we use 
the PSF-subtracted azimuthally averaged radial profile (Fig. 
3) and take the counts in a bin at radial offset of 2.2 arcsec. 
The ratio of this value to the counts in an equivalent bin 
at the position of the star in the unsubtracted data gives an 
estimate of S, the relative brightness of the shell. Note that 
the brightness of the shell at the position of the star is negli- 
gible compared to that of the star itself. We also estimate S 
for the 0.55 WFPC2 observations of Ueta et al. (2000) 
shown in Fig. 3. Since this data is not PSF subtracted. So. 55 
is likely to be an upper limit on the shell brightness in the V 
band. We summarize these quantities, along with the wave- 
length variation, in Table 1. As can be seen, there is re- 
markably little variation in the relative shell brightness, and 
hence the scattering efficiency, between J and K. The ra- 
tio Sj/Sk = 1.32 is much less than the A"'' dependence 
expected for scattering from small (Rayleigh) particles and 
gives a colour excess for the shell of E{J — K) = —0.3. This 
'grey' scattering law was also noted by Kastner & Weintraub 



^ courtesy of the NCSA Astronomy Digital Image Library 
(ADIL) 



0.55 i^m 



Sx X 10^ 1.49 ± 0.09 
Sx/Sk < 3.15 ±0.36 



J 



K 



34.2 ±0.9 32.5 ±1.3 
0.63 ±0.06 0.47 ±0.05 
1.32 ±0.19 1.00 



Table 1. Shell parameters, derived from the PSF subtracted ob- 
servations, in each waveband. P\ is the per cent polarization at 
a radial offset of 1.9 arcsec from the centre of the shell. S\ is a 
measure of the shell brightness relative to the star. The errors on 
Sj and Sk assume a 10 per cent error in the PSF subtraction. 
The error on So. 55 (from which a PSF has not been subtracted) 
is based on photon noise. 



(1995) who suggested that the grains must be large and/or 
have unusual scattering properties. 



2.3.2 Extinction 

An estimate of the extinction to the star can be made based 
on the expected colour for a star of spectral type G5 la and 
the observed fiuxes (HKV89), however there is disagreement 
over the exact stellar temperature. Hawkins et al. (1995) 
found E{B — V) — 0.6, based on a photospheric temperature 
of Tcfi — 5100 K which, assuming a standard interstellar red- 
dening law, leads to Av — 1.8 mag. Reddy & Hrivnak (1999) 
derive a higher temperature of Teb = 6750 K based on an 
analysis of the optical spectrum. The implied bluer intrinsic 
colour leads to a higher estimate of = 4.0 mag. However, 
in a further spectral analysis, Thevenin et al. (2000) dispute 
this higher temperature and derive Tcn = 5660 ±70 K. This 
lower temperature gives a colour excess of E{B — V) = 0.9 
and, again assuming a standard interstellar reddening law, 
an extinction of A^ = 2.9 mag. At a distance of 6 kpc, a sig- 
nificant fraction of this extinction must be interstellar, and 
HKV89 quote 2.0 mag for the interstellar component "as- 
suming a distance of a few kpc". We adopt the constraint 
that the intrinsic extinction, caused by the circumstellar 
shell, should be less than 2 magnitudes in the V band. 



2.3.3 Chemistry 

HD 179821 is an O-rich star with 0/C = 2.6 (Reddy & 
Hrivnak 1999). The envelope chemistry is also O-rich with 
OH masers detected with the same velocity range as the CO 
emission (Likkel et al. 1991, 1989) and lying in a shell of sim- 
ilar extent to the dust shell (Claussen 1993). A weak feature 
at 10 fj,m in the IRAS LRS spectrum (HKV89) indicates the 
presence of silicate grains. More recent ISO results suggest 
that the dust is composed mainly of a mixture of amorphous 
and crystalline silicates (Waters et al. 1996). 



3 SCATTERING MODELS 
3.1 Modelling Procedure 

Details of the Monte Carlo modelling procedure are given in 
Appendix ^. To compare the model results to the observa- 
tions, the model images are convolved with a Gaussian pro- 
file of FWHM 0.9 arcsec, simulating the seeing at the time of 
the observations. They are then processed in the same way 
as the observations, by forming azimuthally averaged radial 
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profiles. This also takes advantage of the symmetry in the 
model to increase the signal to noise of the computations 
for a given number of model photons. The model profiles 
are normalized to the observations using the observed po- 
larized flux peak at 1.9 arcsec, since this is not subject to 
PSF subtraction errors. We calculate the ratio of the shell 
brightness to that of the star for each model and normalize 
this to the observations (5a in Table 1) by varying the dust 
grain number density, A^in. This provides a determination of 
A^in independent of flux calibration, foreground extinction 
or knowledge of the spectrum of the stax. By assuming a 
distance, D, to the object, the physical parameters of the 
model such as optical depth, r, and dust mass, Ma, shown 
in Table 2 are obtained. 



3.2 Spherically symmetric models 

To simulate a spherically symmetric shell we use a dust den- 
sity of the form 



Total Flux 



N{r) = iVi„(r/n„)" 



(2) 



with a = 2, corresponding to a constant mass-loss rate 
throughout the history of the shell. In view of the 0-rich 
chemistry we initially assume a population of bare silicate 
grains. We treat as free parameters the inner radius of the 
shell, Tin, the dust grain number density at this radius, A'in 
and the grain size distribution parameter f3 (see Appendix 
A). In order to generate the observed degrees of polarization 
of 30 to 40 per cent after scattering in a realistic geome- 
try such as a shell, the grains themselves must be able to 
produce 70 to 80 per cent polarization. This requires small 
grains so we initially set P — 6.0 (ASl, Table 2). 

A scattering model using the spherical shell geometry 
and the ASl grain model is shown in Fig. 4. The ring struc- 
ture, seen in polarized flux in the observations (Fig. 1), is 
clearly seen in the model polarized flux image, but not in 
the total flux image. This is because, for lines of sight (LOS) 
that pass through the inner hollow region of the shell, scat- 
tering occurs predominantly from dust in front of the star 
and, therefore, through forward-throwing angles. For these 
angles, there is a large contribution to total flux but little 
to the polarized flux (which is zero for a LOS through the 
star). The effect is to 'fill in' the central region of the shell 
in total flux but to leave a 'hole' in polarized flux. For a 
spherically symmetric detached shell, the detached nature 
of the shell, with a clear boundary at the inner radius rin, 
will only be evident observationally (in scattered light) in a 
polarized flux image. 

The model profiles are shown in comparison with the 
observed profiles in Fig. 5. The observed profiles are plotted 
using the same symbols as in Figs. 2 and 3 but with the 
error bars omitted for clarity. Model calculations at 1.2 fjjn 
and 2.2 /xm are compared with J and K band observations, 
respectively. The characteristics of the model are summa- 
rized in Table 2. We find that the SI model simulates the 
observations quite well, suggesting that a spherical shell is 
a good approximation to the actual geometry of the CSE. 
The shape of the polarized flux profile and the position of its 
peak depend on the scattering asymmetry parameter, g (see 
Appendix B), and the radius of the inner edge of the shell, 
rin. For the SI model, g = 0.34 at 1.2 iim (Table B2). If g 
becomes too large (the grains are too forward-throwing in 
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Figure 4. A spherical shell scattering model (SI) showing the 
distribution of total and polarized flux at 1.2 /im. The images 
have not been smoothed to match the observational data and so 
retain a spatial resolution of 0.18 arcsec. The direct light from 

the star occupies a single pixel at the centre of the image and has 
a total flux count of 764,000 (in arbitrary units). The images are 
scaled between (white) and 50 (black) in the same (arbitrary) 
units. 



their scattering) then the polarized flux peak moves closer 
to the star, as the dust in front of the star begins to domi- 
nate the scattering. We find that, in order to fit the observed 
polarized flux profile, the dust around HD 179821 must have 
g < 0.6. 

To produce the high degrees of observed polarization we 
have had to use a size distribution with a steep power law 
(/? = 6) so that a large proportion of small (highly polariz- 
ing) grains are included in the distribution. This produces 
far too much wavelength dependence in the model and con- 
sequently the shell is much too faint relative to the star at 
the longer wavelength {S1.2/S2.2 = 5.8 in Table 2). This 
suggests that we must construct a model which can produce 
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Table 2. A summary of parameters for each model. Numbers in parentheses are powers of 10. The ASl dust models 
are 'astronomical silicate' with amin = 0-05 /im, Umax = 2.0 fim and /? = 6.0. The AS2 dust model has /3 = 5.5. a, 
e, 7 are parameters of the dust distribution (Section 3.3), rin is the inner radius of the dust shell in m and Ni^ the 
grain number density at that point in m~^. 81,2/82.2 is the ratio of the relative shell brightness at A = 1.2 /im and 

A = 2.2 fim and a measure of the colour of the shell, ri.2 and T2.2 are the optical depths through the shell to the 
star. For the equatorially enhanced models values along both pole and equator (pole/equator) are given. is the 
mass of dust in the shell in Mq. The distance, D, is assumed to be 6 kpc. 



Model 


Dust 


a 


c 


7 


nn 




S1.2/S2.2 


T1.2 


T2.2 


Md 


SI 


ASl 


-2.0 


0.0 


0.0 


1.44(15) 


2.70(-l) 


5.58 


2.57(-l) 


7.70(-3) 


8.20(-2) 


Dl 


AS2 


-2.0 


0.0 


0.0 


1.44(15) 


1.40(-1) 


4.26 


2.27(-l) 


6.68(-2) 


5.10(-2) 


D2 


AS2 


-2.0 


1.0 


10.0 


1.44(15) 


1.14(-1) 


4.22 


1.85(-l)/3.69(-l) 


5.44(-2)/1.09(-l) 


5.70(-2) 


D3 


AS2 


-2.0 


2.0 


3.0 


1.44(15) 


7.50(-2) 


4.10 


1.22(-l)/3.64(-l) 


3.58(-2)/1.07(-l) 


5.90(-2) 


D4 


AS2 


-2.0 


9.0 


3.0 


1.44(15) 


2.90(-2) 


3.93 


4.70(-2)/4.70(-l) 


1.38(-2)/1.38(-l) 


6.60(-2) 


D5 


AS2 


-2.0 


9.0 


10.0 


1.44(15) 


4.60(-2) 


4.09 


7.45(-2)/7.44(-l) 


2.19(-2)/2.19(-l) 


7.20(-2) 



high polarization from larger (greyer) particles. One possi- 
bihty is to move to an axisymmetric dust distribution. 

3.3 Equatorially enhanced models 

The spherically symmetric geometry results in long lines of 
sight through the shell which sample a wide range of scat- 
tering angles. Since maximum polarization is only produced 
for scattering angles axound 90° the geometry acts to reduce 
the maximum polarization of the ASl dust model from ~ 80 
per cent to ~ 40 per cent. By concentrating the dust in the 
plane of the sky, scattering angles around 90° are given more 
weight and the polarization can be increased. This can be 
achieved using a model with an equatorial density enhance- 
ment which is then viewed pole-on. To simulate this we use 
a dust density of the form 

iV(r) = iVi„(r/nn)-"(l-hesin'^6>) (3) 

which is similar to that used to generate axisymmetry in 
models of planetary nebuleie (e.g. Kahn & West 1985; Kwok 
2000) . The radial density fall off is now modified to create a 
density enhancement of 1 + e between the equator and pole. 
The parameter 7 controls the way density increases with 
polax angle, 9, with higher values of 7 concentrating more 
material in the equatorial region. 

In Fig. 6 we show plots of the radial dependence of 
surface brightness and polarization at 1.2 fim for five 'disc' 
models with various values of e and 7. These models, Dl to 
D5, simulate increasing concentrations of dust in the equato- 
rial plane, with Dl being a spherically symmetric model and 
D5 being the most flattened with a density contrast (equa- 
tor to pole) of 10:1 and an angular dependence of sin^° ^. 
All are viewed pole-on and have a = 2 with the inner ra- 
dius, Hn, fixed at 1.44 x 10^® m (D = 6 kpc). A sihcatc dust 
model, AS2, is used which difi'ers from the ASl model used 
in Section 3.2 in that it has a slightly shallower power law 
for the size distribution with /3 = 5.5, thus including larger 
grains. 

Only the D5 model with a 10:1 density contrast (e = 10, 
7 = 10) approaches the observed polarization. However, in- 
creasing the density contrast this much also acts to flatten 
the polarization profile so that it no longer matches the ob- 
served rise in polarization with radial offset. Increasing the 
density contrast between equator and poles also causes the 
intensity profile to turn over at small radial offsets until, for 
the D5 model, there is a peak in intensity similar to that 
observed in polarized flux. In other words, there is now a 



central 'hole' in intensity as well as polarised flux (Fig. 1). 
There is no evidence for such a hole or inner shell boundary 
in the V-band HST images (Ueta et al. 2000), although run- 
ning the D5 simulation at 0.55 /im also produces a central 
hole. This indicates that any density enhancement must be 
slight (Dl and D2 models). 

It is also apparent from Table 2 that concentrating the 
dust in the equatorial plane does nothing to solve the colour 
problem - the shell is still too blue. The S1.2/S2.2 index does 
decrease slightly for the more axisymmetric models but this 
is due to an increase in extinction for paths in the equatorial 
plane. Changes to the geometry alone cannot increase the 
polarization produced by the dust models based on astro- 
nomical silicate so as to simultaneously match the observed 
polarization and still retain colour neutral scattering at J 
and K. 



4 DISCUSSION 

The results of the previous section indicate that the radial 
distribution of polarized flux, surface brightness and per cent 
polarization in the near-IR can be modelled using a spheri- 
cal distribution of dust with a sharp inner boundary and an 
density fall off. It has not been possible to simultane- 
ously fit the degree of polarization and the low colour excess 
{E{J — K) = —0.3) since all of the models considered pro- 
duce scattering which is far too blue. The shell could easily 
be reddened by increasing the optical depth, however, this 
would also reduce the polarization. In addition, examination 
of the V band HST image of Ueta et al. (2000) suggests that 
the shell is optically thin in the V band, since the surface 
brightness profile (Fig. 3) is very similar to that of our J 
band data. 

So far our dust model has been based on 'astronomi- 
cal silicate' (Draine 1995) which has been widely used to 
model reflection nebulae in stax formation regions. How- 
ever, silicates only produce high degrees of polarization for 
small grains (close to the Rayleigh limit) which inevitably 
exhibit strong colour effects in the scattered intensities, pro- 
ducing predominantly 'blue' reflection nebulae. In Appendix 
B we tabulate the scattering properties of various size dis- 
tributions of 'astronomical silicate' and other materials. For 
an optically thin nebula (where absorption is not impor- 
tant), the ratio of scattering cross sections at 1.2 and 2.2 /im 
gives a reliable indication of the nebula colour. In order to 
produce a nebula with a ratio of scattering cross sections 
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Figure 5. A comparison of azimuthally averaged radia 
through the spherical slicll model SI with those through 
(Figs. 2 and 3). The error bars, and data points in th 
region, have been omitted for clarity. The model profiles a 
to the observations by normalizing the model polarized flux (solid 
line) to the observed polarized flux (x) at a radius of 1.9 arcsec. 
This also determines the scaling on the total flux curves (dashed 
lines). Per cent polarization is shown as a dotted line. 
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Figure 6. Azimuthally averaged radial profiles through shell 
models with equatorially enhanced dust distributions, compared 
with the data (Figs. 3). The error bars on the data points have 
been omitted for clarity. Profiles of both total flux (upper panel) 
and per cent polarization (lower panel) are shown. The parame- 
ters for the models are listed in Table 2. 



Csca(1.2/im)/Csca(2.2/im) < 2.0, a power law size distri- 
bution of astronomical silicates with Omin = 0.05 /um and 
Omax = 2.0 /im must have an index /3 < 4.0. Such a distri- 
bution produces very little polarization and cannot fit our 
observations. The same is true of the other forms of silicate 
listed (olivine, pyroxene) which have similar refractive in- 
dices in the near-IR. As mentioned in Section 2.2, the PSF 
subtraction is only approximate, so that it is possible we 
have over-estimated the intrinsic shell polarization, in which 
case some of the more highly polarizing grain models (Ta- 
ble B2) may account for the observations. However, even if 



we take the raw (before PSF subtraction) polarization of 
20 per cent at J and K as the intrinsic polarization of the 
shell, a bare silicate dust model cannot produce the required 
colours. We therefore rule out scattering from spherical sili- 
cate particles as the dominant polarizing mechanism. 

4.1 Dust composition 

Some of the emission features seen towards HD 179821 may 
be due to metal oxides in the shell, such as magnesium- 
iron-oxides (Ueta et al. 2001). There is evidence that these 
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materials may contribute to the observed structure in the 
15-24 jj,m spectral region (Henning et al. 1995, Waters et 
al. 1996). We have calculated the scattering properties for 
Feo.4Mgo.6O (Henning et al. 1995) and find that this mate- 
rial produces more polarization than silicates, for an equiv- 
alent size distribution. A scattering cross section ratio of 
Ri.2/2.2 < 2.0 can be produced for a power law size index of 
/3 < 4.5 (Table B2) which results in maximum polarization 
of < 60 per cent at 1.2 nm and < 50 per cent at 2.2 /itm. This 
is still not sufficient to fit the observed shell polarization, al- 
though it could account for the ~ 20 per cent polarization 
observed in the raw (non-PSF corrected) data. The mate- 
rial docs produce around 10 per cent more polarization at 
1.2 nm than at 2.2 /irn, though, which is contrary to our 
J and K band observations. In addition, even if metal ox- 
ides are present in the shell, it is not clear that they would 
be majority contributors to the scattered light or that they 
would be distributed so evenly. 

The ISO spectra of HD 179821 also show evidence for 
crystalline ice in the form of a 43 jj,m 'bump', along with 
various features attributable to forms of crystalline silicates 
(Waters et al. 1996) . We have calculated the scattering prop- 
erties for ice coated silicates with a coating thickness equal 
to the grain radius (Table B2). Ices, especially water ice with 
a low refractive index, are highly polarizing so that a coat- 
ing of ice can increase the polarization of a silicate grain 
core significantly (by an amount increeising with the coat- 
ing thickness). However, the scattering properties of the ice 
are more wavelength sensitive than silicates, so that there 
is a trade off between increasing the maximum polarization 
of the dust by applying an ice coating and maintaining a 
grey scattering efficiency. Again we have not been able to 
produce the required high polarization and grey scattering 
using ice coated silicates. 

4.2 Particle size distribution 

The power law size distribution used in the previous sec- 
tion has upper and lower cut-offs on the grain size so as to 
keep the number of grains finite. We have adopted the values 

cimin=0.05 /im and amax=2.0 fiia, rather than treating them 
as free model parameters, and so we should examine these 
assumed values. For any size distribution, the particle sizes 
that contribute most to scattering will be those that max;i- 
mize the weighted scattering cross section, N{a)Csca{a). In 
order to produce the high polarizations that we see in the 
J and K band observations, these particles must be small, 
typically with size parameter a; < 1 (or radius a < 0.2 /xm 
at A = 1.2 j-im), where x = 2na/\. This is illustrated in 
Fig. 7 which shows the polarization produced at A = 1.2 /xm 
as a function of a. Particles larger than x = 1 produce less 
polarisation and of varying sign (especially for the dielectric 
materials such eis silicate) which sum to small values when 
a distribution of sizes is taken. Since the scattering cross- 
section increases with particle radius, a, as up to the sixth 
power, a steep power law index, /3, is required for the size 
distribution to ensure that the smaller (polarizing) parti- 
cles arc present in sufficient numbers to dominate the scat- 
tered light. This in turn means that the scattering properties 
are rather insensitive to the upper size cut-off, Omax, since 
there are too few particles of this size to contribute. At the 
small size end of the distribution (a; < 1), scattering is ap- 



proaching the Rayleigh limit where the cross-section falls 
off as Csca oc a® and polarization is a maximum. As long 
as flmin IS small enough to include Rayleigh particles, then 
reducing it further has very little effect on the polarization 
or colour of the scattered light. For values of /3 < 6, the 
smallest particles do not contribute to the scattered flux. 
Fig. 7 shows that it is not possible to generate the ob- 
served high polarization and neutral nebula colour by ad- 
justments to the power law size distribution, at least for the 
materials we have considered. The range of particle sizes re- 
sponsible for producing high polarization and neutral colour 
(low Csca(1.2/nm)/Csca(2.2//m)) are mutually exclusive. This 
would also be the case for other size distribution functions, 
such as a modified power law (Jura 1994) or exponential 
distribution (Jura 1975). 

4.3 Non-spherical or composite particles 

It is well known that much of the dust in the Galaxy is 
composed of non-spherical grains which are responsible for 
phenomena such as the interstellar linear and circular po- 
larization (Whittet 1995; Martin 1972). Mid-IR polarization 
observed in star-forming regions (e.g. OMCl, Aitken et al. 
1997) shows that aligned, non-spherical particles are present 
in circumstellar environments also. Given the wide-spread 
distribution of non-spherical grains it seems reasonable to 
assume that at least a proportion of the dust at its point of 
origin, in the circumstellar envelopes of evolved stars, will 
also be non-spherical. 

The scattering properties of highly spheroidal particles 
differ significantly from spheres and can show a combina- 
tion of the characteristics of larger and smaller particles. 
Zakharova & Mishchenko (2000) find that for spheroids with 
aspect ratios of 20:1, the polarization reaches 100 per cent 
at a scattering angle of 90 ° , typical of Rayleigh scattering, 
even for particles with size parameters x = 3.5. However, 
the scattering efficiency as a function of size parameter is 
more typical of that of larger spheres, showing a weak de- 
pendence on X. It therefore seems plausible that the high 
polarizations and weak dependence of Csca on grain size, re- 
quired to explain our observations, may be satisfied if the 
scatterers are spheroidal with large aspect ratios (needle- 
like or plate-like). An alternative to highly elongated parti- 
cles could be provided by composite or aggregate particles 
where the overall size of the particle is comparable to the 
wavelength, but many of the scattering properties (such as 
the polarization) are characteristic of the component sizes. 
Such behaviour has been seen in low density aggregates of 
spheres (West 1991). 

4.4 Dust mass, extinction and the distance to HD 
179821 

The dust optical depths shown in Table 2 have been obtained 
by fitting to the scattered surface brightness of the shell rel- 
ative to the star, and are intrinsic to the object and not 
dependent on the assumed distance. If the object is closer 
than 6 kpc then path lengths through the shell would de- 
crease and the dust density would have to be increased to 
maintain the observed shell brightness (and hence optical 
depth). So our model fits do not directly constrain the dis- 
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Figure 7. Linear polarization at A = 1.2 /im (top panel) and the ratio of scattering cross sections at A = 1.2 and A = 2.2 fj,m 
(bottom panel) as a function of grain radius, a, for various materials. 



tance estimates. Using the silicate dust model and spheri- 
cal (SI) geometry we calculate an extinction to the star at 
0.55 jj,m (roughly V band) of ro.s.s = 1.67. This is consistent 
with estimates of the extinction to the star in the V band 
of between 1 and 2 mag (Section 2.3.3). 

For the SI model we calculate a dust mass of Ma = 
0.08 Mq, assuming an outer shell radius of 9 arcsec (Kast- 
ncr & Wointraub 1995), a distance of 6 kpc and a density 
of 3.0 X 10^ kg m"'' for silicate dust. This is of the same 
order as estimates of the dust mass obtained from mid-IR 
modoHing of 0.03 M© (Jura & Werner 1999) and 0.04 Mq 
(Hawkins et al. 1995). These authors assume values for the 
gas-to-dust mass ratio of between 100 and 160 which gives a 
total (gas-|-dust) mass in our model of between 8 and 13 M©. 
Our estimate of Ma is sensitive to the dust model, in partic- 
ular the scattering cross-section, and could easily be smaller, 
in agreement with the estimates from mid-lR emission. For 
example, increasing the scattering cross-section by using 
slightly larger grains (Dl in Table 2) gives Ma = 0.05 M©, 
although this no longer fits the observed polarizations. 



5 CONCLUSIONS 

The circumstellar envelope of the evolved star HD 179821 
is seen by scattered Ught in the near-IR and appears as a 



detached shell in linearly polarized light. We have used ob- 
servations of a PSF calibration star to subtract the atmo- 
spherically scattered PSF of HD 179821 in order to obtain 
an image of the shell and hence calculate its intrinsic linear 
polarization. We find high degrees of polarization, between 
30 and 40 per cent, in both J and K bands, with evidence 
for a gradual increase in polarization with radial offset from 
the star. The colours are remarkably neutral for scattering, 
with a colour excess in the shell of E{J — K) = —0.3. 

We use an axisymmetric scattering code to model the 
J and K band observations and find that the distribution 
of polarized and total flux can be produced by a spherically 
symmetric shell of dust with an r^^ density distribution. 
The model shows that, when scattering from such a spheri- 
cally symmetric shell, the detached nature of the shell is only 
apparent in polarized flux, with the total flux image showing 
no evidence for an inner boundary. By fitting the models to 
the observed brightness of the shell relative to the star, we 
obtain an optical depth to the star at 1.2 /im of ti.2 = 0.26 
which is consistent with estimates of the V band extinction. 
Assuming a distance to HD 179821 of 6 kpc, we derive an 
inner shell radius of ri„ = 1.44 x 10^^ m, a dust number 
density of Nin = 2.70 x 10~^ m~^ at rfn and a total mass 
of dust in the shell of Ma = 0.08 M©, assuming a power- law 
size distribution of silicate grains. These estimates are sen- 
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sitive to uncertainties in the dust model but the dust mass 
is of the same order as that derived from mid-IR imaging of 
the warm dust. 

We have experimented with equatorially enhanced den- 
sity distributions, viewed polc-on, but conclude that if the 
shell does have an increased concentration of dust in the 
equatorial plane then the contrast between equator and 
poles cannot be much more than 2:1, otherwise the radial 
dependence of linear polarization becomes too flat to fit the 
observations. 

Although the spherical shell model fits the observed dis- 
tributions of total and polarized flux quite well, we have not 
been able to simultaneously account for the high degrees of 
linear polarization and the unusually 'grey' scattering at J 
and K. We attribute this to inadequacies in the dust model 
and have explored various materials such as silicates, car- 
bons, metal oxide and also ice-coated particles. If the high 
intrinsic polarizations implied by our PSF subtraction are 
correct, then none of the materials has reproduced the ob- 
served colours and polarizations and a distribution of homo- 
geneous spherical particles, as assumed here, is unlikely to do 
so. Instead, we suggest that the dust surrounding HD 179821 
may consist of either highly elongated noodle- or plato-likc 
particles or of composite grains with many smaller inclu- 
sions. 

It woxild bo very useful to obtain higher spatial res- 
olution imaging polarimctry observations (such as corono- 
graphic 8-rn data) and/or spcctropolarimctry over the op- 
tical and near-IR, to ascertain the polarization in the shell 
more accurately and to determine the scattering behaviour 
of the dust over a wider wavelength range. 
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APPENDIX A: DETAILS OF THE MONTE 
CARLO SCATTERING MODEL 

The scattering models are computed using a modified ver- 
sion of the Monte Carlo scattering code of Menard (1989). 
A given number of photons (usually lO'^ for the optically 
thin models) are followed as they leave the star and propa- 
gate through the dusty envelope. At each scattering event, 
the new photon path and Stokes intensities are calculated. 
The density structure of the envelope can be spherically 
symmetric or axisymmetric (2D) and is specified as a ta- 
ble of grain number density at each point (r, 9) where r 
is the radial distance from the star and 9 is the polar an- 
gle (9 = for a pole-on view and 9 = 90 for an equatorial 
view). Photons emerging from the nebula are collected in 
bins according to their polar angle to form a 2-D scattered 
light image of the nebula for each viewing inclination angle. 
In practice these bins are averaged over inclination angle, 
to improve signal to noise, so that the bin approximating a 
pole-on view (bin 1) actually encompasses viewing angles in 
the range 61 = 10°. Since the observations of HD 179821 
suggest that the nebula geometry is either spherically sym- 
metric, or axisymmetric but viewed pole-on, this bin is used 
throughout. 

The dust grain model is specified by the refractive index 



© 0000 RAS, MNRAS 000, 000-000 



The dust shell around HD 179821 



11 



of the dust, m, and the grain size distribution. Wc assume 
a truncated power law distribution of grain sizes such that 
the number of grains of size a per unit volume is given by 
n(a) = Aa~'^ where amin < a < amax and ^4 is a constant 
defined such that 



and TV is the total number of grains per unit volume, which is 
given at each {r,6) position by the density table. The grains 
are assumed to be spherical and the scattering calculations 
use Mie theory. 

The model geometry suggested by the observations is 
a hollow shell, with inner radius rin and outer radius rout. 
The number density of dust grains at the inner radius, rin, 
is A'in, and outside the shell (r < rin and r > rout) it is 
zero. Within the shell the grain number density is given 
by the density table. We set the outer shell boundary to 
5.4 X 10'' AU, corresponding to 9 arcsec at 6 kpc, as indicated 
by the extent of the scattered light detected by Kastner & 
Weintraub (1995). The model images at each viewing angle 
are computed on a 2D grid of 99 x 99 pixels, with the star 
at the central pixel. For modelling purposes we have ignored 
the observed offset of 0.3 arcsec between the star and the 
centre of the shell. With an outer shell radius of rout = 9 
arcsec then a model pixel is equivalent to 0.18 arcsec. 



APPENDIX B: SCATTERING PROPERTIES OF 
VARIOUS DUST GRAIN MODELS 

We list here the results of scattering calculations for various 
particle size distributions. All calculations, assume Mie scat- 
tering and are based on modified versions of the BHMIE and 
BHCQAT spherical particle algorithms of Bolircn & Huffman 
(1983). Power-law size distributions (Section 3.1) are used 
throughout with Omin = 0.05 /xm and ttmax = 2.0 /xm. Table 
Bl lists optical constants at A = 1.2 iiui and A — 2.2 fim for 
various materials and Table B2 lists their scattering proper- 
ties. 

The albedo, uj is the ratio of scattering and extinction 
cross-sections. The asymmetry parameter, g, is the average 
cosine of the scattering angle and a measure of the forward- 
throwing nature of the scattering. For a particle which 
scatters light equally in all directions, g = 0, whereas for 
forward-throwing scattering, g lies in the range < p < 1. 
More formal definitions may be found in Bohren & Huffman 




(Al) 



(1983). 
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Table Bl. Optical constants for various materials at A = 1.2/ini and A = 2.2/im. 



Material 


m(A = 1.2/im) 


m(A = 2.2/im) 


Comment 


Ref 


Ast Sil 


1.71, 0.03 


1.71, 0.03 


'astronomical silicate' 


1 


Mgo.8Fei.2Si04 


1.790, 0.103 


1.866, 0.052 


Fe-rich olivine glass 


2 


Mgo.6Feo.4Si03 


1.640, 0.004 


1.620, 0.002 


Fe-rich pyroxene glass 


2 


Feo.4Mgo.6O 


1.95, 0.31 


2.15, 0.26 


iron-magnesium oxide 


3 


Amor. Carbon 


2.51, 0.76 


2.79, 0.80 


amorphous carbon 


4 


H2O Ice 


1.298, 6.7(-6) 


1.262, 2.6(-4) 


water ice at 140 K 


5 



References: 1 Draine 1985; 2 Dorschner et al. 1995; 3 Henning et al. 1995; 4 Preibisch et al. 1993; 5 Warren 1984. 



Table B2. Scattering properties for power-law size distributions of spherical particles. /3 is the power-law size 
distribution index, Csc<i(1.2) and Csca (2.2) are the scattering cross-sections at 1.2 and 2.2 //m, 0)1.2 and 0^2.2 the 
albedos, gi,2 and §2.2 the scattering asymmetry parameters, Pm(1.2) and P^(2.2) the maximum lineaj: polarizations 
and Ri.2/2.2 is the ratio Csca(1.2)/Csca(2.2). 



Material 


P 


sea 


1 o^ 
1.2) 


^ sea 


2.2) 


Ast Sil 


6.0 


5.27( 


-16) 


9.071 


-17) 




5.5 


1.02( 


-15) 


2.291 


-16) 




5.0 


2.15( 


-15) 


6.351 


-16) 




4.5 


4.88( 


-15) 


I.9O1 


-15) 




3.5 


3.02( 


-14) 


1.961 


-14) 


TV T ■ /" A 

Mgo.8Fei.2Si04 


6.0 


5.69( 


-16) 


1.221 


-16) 




5.5 


1.06( 


-15) 


2.971 


-16) 




5.0 


2.12( 


-15) 


7.9O1 


-16) 




4.5 


4.53( 


-15) 


2.221 


-15) 




3.5 


2.53( 


-14) 


2.041 


-14) 


Mgo. 6*60.40103 


6.0 


4.66( 


-16) 


7.6I1 


-17) 




5.5 


9.31( 


-16) 


1.981 


-16) 




o.U 




-15) 


0. r 01 


-lb) 




4.5 


4.83( 


-15) 


1.80( 


-15) 




3.5 


3.31( 


-14) 


2.06( 


-14) 


Feo.4Mgo.6O 


6.0 


6.71( 


-16) 


1.45{ 


-16) 




5.5 


1.16( 


-15) 


3.27( 


-16) 




5.0 


2.17( 


-15) 


8.00( 


-16) 




4.5 


4.34( 


-15) 


2.14( 


-15) 




3.5 


2.17( 


-14) 


1.66( 


-14) 


Amo Caj" 


6.0 


1.11( 


-15) 


2.16( 


-16) 




5.5 


1.80( 


-15) 


4.50( 


-16) 




5.0 


3.08( 


-15) 


1.01( 


-15) 




4.5 


5.66( 


-15) 


2.42( 


-15) 




3.5 


2.46( 


-14) 


1.71{ 


-14) 


H2O ice 


6.0 


1.20( 


-16) 


1.56( 


-17) 




5.5 


2.63( 


-16) 


4.30( 


-17) 




5.0 


6.63( 


-16) 


1.37( 


-16) 




4.5 


1.89( 


-15) 


4.80( 


-16) 




3.5 


1.97( 


-14) 


7.11( 


-15) 


Mgo.8Fei.2Si04 H 


h H2O ice 6.0 


1.72( 


-16) 


2.68( 


-17) 




5.5 


3.60( 


-16) 


7.29( 


-17) 




5.0 


8.45( 


-16) 


2.23( 


-16) 




4.5 


2.20( 


-15) 


7.41( 


-16) 




3.5 


1.92( 


-14) 


9.58( 


-15) 



t a negatively polarised backscattering peak is also produced 







91.2 


92.2 


-Pm(1.2) 


-Pm(2.2) 


Ri.2/2.2 


0.62 


0.37 


0.34 


0.32 
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